Optimal entanglement generation in cavity QED with dissipation 



O 

o 

O 

Q 

CO 
CN 

S 1 

cr 



> 
o 

CN 

o 



X 



Jian-Song Zhang 1 ' 2 and Jing-Bo Xu^ 

1 Zhejiang Institute of Modern Physics and Physics Department, 
Zhejiang University, Hangzhou 310027, People's Republic of China 
2 Department of Applied Physics, East China Jiaotong University, 
Nanchang 330013, People's Republic of China 

Abstract 

We investigate a two-level atom coupled to a cavity with a strong classical driving 
field in a dissipative environment and find an analytical expression of the time evo- 
lution density matrix for the system. The analytical density operator is then used to 
study the entanglement between the atom and cavity by considering the competing 
process between the atom-field interactions and the field-environment interactions. 
It is shown that there is an optimal interaction time for generating atom-cavity 
entanglement. 

OCIS codes: 270.5580, 270.5585 



I. INTRODUCTION 

Quantum entanglement is one of the most 
striking features of quantum mechanics [j], 0|. 
It is central to many active research fields, 
such as quantum computation (3|, quantum 
teleportationQ], quantum communication [Ej], 
and quantum key distribution [§[. The gen- 
eration of quantum entangled states has at- 
tracted much attention and many physical 
systems including cavity quantum electrody- 
namics (Cavity QED)@, 1, 1, [To), trapped 
ions[ll||, nuclear magnetic resonance [II and 
quantum dots [l^], have been suggested to 
generate entangled states. 

Solano et 



In Ref. 14 



al. proposed a 
scheme to generate the Schrodinger cat states 
using two-level atoms interacting with a cav- 
ity with the help of a strong classical driv- 
ing field. However, a quantum system is 
unavoidably influenced by its surrounding 

The interaction 



environment 13, 16, 17, 18 



between the quantum system and its environ- 
ment leads to decoherence which is the main 
problem in the generation of quantum entan- 



glement. In[19|, the authors investigated the 



nonclassical properties of a cavity field when 
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no atom crosses it under the influence of dis- 
sipation by making use of phase-space tech- 
niques. 

In the present paper, we consider a quan- 
tum system consisting of a two-level atom 
coupled to a cavity of high quality factor with 
a strong classical driving field in a dissipative 
environment and find an analytical expres- 
sion of the time evolution density matrix for 
the system. The analytical density matrix is 
then used to study the entanglement and the 
purity of the system. In this model, there 
are two competing processes in this system. 
One is the amplification process introduced 
by the atom-field interactions (which leads to 
entanglement). The other is the dissipation 
process due to the field-environment interac- 
tions( which leads to disentanglment). Our 
calculation shows that there is an optimal in- 
teraction time for generating atom-cavity en- 
tanglement when the coupling constants and 
the frequencies of the two-level atom, the cav- 
ity, and the classical driving field are given. 
In section II, we solve the master equation of 
the system by making use of the superoper- 
ator algebraic method and obtain an explicit 
expression of time evolution density matrix 
for the system of a two-level atom coupled to 
a cavity with a strong classical driving field 



2 



in a dissipative environment. In section III, 
we use the concurrence to investigate the en- 
tanglement between the two- level atom and 
the cavity by means of the analytical expres- 
sion of the density matrix for the system. In 
section IV, we calculate the purity of the sys- 
tem by employing the linear entropy. A con- 
clusion is given in section V. 



II. SOLUTION OF AN ATOM IN A 
DECAY CAVITY WITH A STRONG 
CLASSICAL DRIVING FIELD 

Now, we consider the time evolution of a 
two-level atom driven by an external classical 
field in a dissipative cavity. The total Hamil- 
tonian of the system reads 14J 



H 



i ( i \\ 

Y a z + go{cr+a + cr-a 'j 



+A (o-+e 



<j_e 



where u, u and u c are the frequencies of 
the cavity, atom and classical field, respec- 
tively. The operators a z and o± are defined 
by o z = |0> <0| - |1)(1|, a + = |0)(1|, and 
er_ = |1)(0| where |0) and |1) are the excited 
and ground states of the atom, respectively. 
Here, a and a) are the annihilation and cre- 
ation operators of the cavity; go and Ao are 
the coupling constants of the interactions of 
the atom with the cavity and with the clas- 
sical driving field, respectively. Note that we 
have set fi = 1 throughout this paper. For 
the sake of simplicity, we consider the res- 
onant interaction between one mode of the 
cavity and the two-level atom strongly driven 
by the classical external field. The schematic 
picture of the present system may be modeled 
as shown in Fig.l. 

In the strong driving regime Ao go, 
we can realize a rotating-wave approximation 
and eliminate from Eq.(l) the terms that os- 
cillate rapidly. The Hmailtonian of the sys- 
tem in the interaction picture is 14 




cavity 



Classical driving field 

FIG. 1: The schematic picture of a two- level 
atom driven by an external classical field in a 
dissipative cavity. 

with g = go/2. We assume that there is a 
reservoir coupled to the electromagnetic field 
in the cavity. Thus, the master equation that 



governs the dynamics of the system is 15 



| = -Hv, p]+ v P 



-i[V, p] + k{2apoJ — cvap — pa^a), (3) 



where k is the decay constant and the super- 
operator T> = k(2a-a^ —a)a — a^a) represents 
the losses in the cavity. 

Now, we show how to obtain an analytical 
solution of the master equation (3) with the 
help of superoperator algebraic method [2l|, 



251 ] . The density matrix for the system can 



be represented as follows 

p(t) = Poo (t)® |0)(0|+p n (t)®|l)(l| 

+P0i(*)® |0)(l|+p 10 (f)®|l)(0|, (4) 

where pi/s are defined as pij = (i\p\j), i,j = 
0,1. A straightforward calculation shows 
that 

Poo = -ig{{a + a f )pio - Poi(a + a r )} + Vp 00 , 
Poi = -ig{{a + a f )pn - PooO + a r )} + Vp 01 , 
Pio = ~ig{(a + a 1 ) poo - p u (a + a f )} + Vp w , 
Pn = -ig{(a + a f )poi - PioO + a 1 ")} + "Ppn- 



V = g(a + a^(a + + a_), 



(2) 



(5) 
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It can be proved that the above equations can use of the superoperator algebraic method 
be recast as 
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v±(t) 
i±{t) 

with 



C>'±rfJt)* 



(6) 
(7) 



C± = ±ig(a ■ +<r ■ — ■ a — -a') + T>, 

C'± = ±ig(a ■ +a t ■ + ■ a + -a 1 ") + V, 

V± = (Poo + Pn) =F (poi + Pio), 

V± = (Poo - Pn) T (Pio - Poi)- (8) 

Here, the superoperators a-, -a, a/- and -a) 
represent the action of creation and annihila- 
tion operators on an operator 



(a-)o = ao, (-a)o = oa, 



(9) 



Combing the communication relation 
[a, a'] = 1 with the above equation, we 
derive the following communication relations 
between the superoperators 



[(a-),(at-)] = l, [(-a), (-at)] 



-1, 



(10) 



while all other communication relations are 
zero. It is easy to check that the superoper- 
ators M. = a ■ a) , 1Z = a) a- and £ = -a) a 
satisfy the relations 

[K,M} = [£,M] = -M,[K,C}=0. (11) 

We assume the field is initially in coherent 
state \a) and the atom is in the exited state 
|0). Thus the initial state of the atom-field 
system is 



|V(0)> = |a>®|0). 



:i2) 



Then the elements of the density matrix 
are initially poo(0) = |c*)(a| and pn(0) = 
Poi(0) = Pio(0) = which leads to rj±(0) = 
rf ± {0) = \a)(a\. 

Using the above results we find the solu- 
tion of the master equations (6) by making 



v±(t) 



c±t 



»7±(0) 



e ±i9tyi+ktXl {\a){a\) 

e ^ le E£(l-e- W(|a)(a|) 
e {e 2kt -l)M e -ktn e -ktC 



|a±(t))(a±(t)|, 



(13) 



where 

\a±(t)) 
Xi 



\ae 



e' kt )), 



2M-TZ-C, 

a ■ +a) ■ — ■ a — -a) . 



(14) 



Similarly, the solutions of the master equa- 
tions (7) can also be derived as 



±igt{X2+y 2 )+ktX 1 



e J e 



a)(a\ 
a){a\ 



f±(t)\a ± (t)){a T (t)l 



(15) 



where 



h 

h 

f± 



2(a • + ■ a) 

-v 



a ■ +a) • — ■ a' + -a, 



kl 



1 + kt), 



k 2 v 
|d-e-), 

J l +2\f 2 \ 2 ±h{a+a*)(2-e- kt ) 



(16) 



Combing the solutions of the Eqs.(6) and 
(7), we obtain the explicit expression of the 
density matrix of the system at time t 

p(t) = p o(t)® |0)(0|+pn(t)(8)|l)(l| 

+P01W® |0)(l| + pi (t)<8)|l)(0|,(17) 
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with 

PooO) 

Pn{t) 
Pw{t) 



-[\a + (t))(a+(t)\ + \a-(t))(a4t)\ 

+f + (t)\a + {t))(a-(t)\ 
+/_(t)|a_(t))(a 



-[|a + (t))(a + (t)| + |a_(t))(a_(t)| 
-/ + (f)|a+(f)>(a_(t)| 
-/_(t)|a_(f))<a + (f)IL 
~[\a+(t))(a+(t)\-\a4t))(a-(t)\ 

-f + (t)\a + (t))(a-(t)\ 
+/_(t)|a_(t))(a + (t)|], 

p!i(*)- (18) 



Thus, we find an analytical solution of 
the master equation with the initial state 
\i/j(0)) = \oc) (g> |0). In the following sections, 
the explicit expression of density matrix is 
used to investigate the entanglement and the 
purity of the system. 



III. ENTANGLEMENT BETWEEN A 
TWO-LEVEL ATOM AND DECAY 
CAVITY 

It is well known that entangled states 
are the basic resource of quantum informa- 
tion processing, such as quantum communi- 
cation and quantum teleportation. Cavity 
QED is a useful tool to generate entangled 
states. However, an entangled state will be- 
come mixed and/or less entangled due to the 
decay of cavities. Here, we confine our con- 
sideration in the question whether we can 
choose an optimal combination of parameters 
g, k, and t to maximize the entanglement of 
the atom-field system. In order to study the 
entanglement of above system described by 
density matrix p, we adopt the measure con- 
currence which is defined as 1221 



C = max {0, Ai — A 2 — A 3 — A 4 }, 



(19) 



where the A, (i= 1,2, 3, 4) are the square roots 
of the eigenvalues in decreasing order of the 



ay ) and 



(T y is the Pauli Y matrix, i.e., a y 



magnitude of the "spin-flipped" density ma- 
trix operator R = p{a y ® cr y )p*(a y (g> 

'0 -i 
i 

In general, the model consisting of the 
two- level atom and the field is a 2 x oo system. 
However, for the initial state ^(O)) = \a) ® 
|0), the atom- field system can be mapped 
onto a 2 x 2 system as one can see from the 
density matrix. To this end we introduce two 
orthonormal vectors | |) and | j) which are 
defined by 



T) 
I) 



| «+(*)), 



Vl-|A(*)| : 



■(!«-(*)> 



(20) 

A(t)K(t)», 



with A(t) = (a + (t)\ a -(t)) ■ Therefor, the 
states \a + (t)) and |a_(£)) can be represented 
in terms of | |) and | |) 



M*)> 

|a_(t)> 



IT), 

A(t)it) + v^ 



|A(t)| 2 |!>. (21) 



Inserting Eq.(21) into Eq.(18), one can 
rewrite the state of the atom-field system 
at time t, i.e., the state of system can be 
mapped onto a system composed by two two- 
dimensional subsystems. Note that the den- 
sity matrix of Eq.(17) is now an effective 
two-qubit system whose entanglement can be 
evaluated by the measure concurrence. With 
the help of Eqs.(17), (18), (19), and (21), one 
can obtain the concurrence the density ma- 
trix of Eq.(17) which is too long to present 
here. In Fig. 2, we plot the entanglement dy- 
namics of the quantum system as a function 
of time t. It is worth noting that we have used 
the dimensionless quantities throughout this 
paper. 

We now turn to show that our model can 
be used to generate the Schrodinger cat state 
in the case of a = k = 0. A simple calculation 
shows that 



a±{t) = ±igt,f 1 (t) = - 

Mt) = igt,f ± = l, 

V±(t) = I ±igt)(±igt\. 

rj' ± (t) = \ ±igt)(Tigt\. 



-2gH\ 



(22) 



5 



Here, we have used the fact a = k = 0. Thus 
the atom-field state at time t is 

\m) = - igt) + \igt)) ® \o) 



+{\-igt)-\igt))®\l)\ 

A=(\-igt)\+)-\igt)\-)), 



with 



1 



l±> = ^ID±|0>. 



(23) 



(24) 



The entangled state \ift(t)) is the Schrodinger 
cat state. 

In the case of k > 0, entangled states will 
become mixed and/or less entangled due to 
the decay of cavities. Fortunately, one can 
find that there is an optimal time t opt to max- 
imize the entanglement if the coupling con- 
stant g and the decay rate of the cavity k are 
given. In Fig.2, we plot the concurrence of 
the atom-field system as a function of time t 
and the coupling constant g for k = (Upper 
panel) and k = 0.05 (Lower panel). It is not 
difficult to see that in the case of k = the 
concurrence first increases with time t and 
then reaches a plateau. However, when k > 
there is no plateau for the concurrence. It 
first increases with time t, and reaches the 
maximal value, then decreases with time t. 
Finally, the atom and the field is disentan- 
gled. In order to see this more clearly, we 
plot the concurrence as a function of time for 
different values of decay constant k in Fig.3. 
One can easily find that if the coupling con- 
stant g and the decay rate of the cavity k are 
known one can find an optimal time t op t to 
maximize the concurrence. For instance, the 
optimal time t opt are 0.83 and 0.61 in the case 
of g — 1, k — 0.1 and g — 1, k — 1, respec- 
tively. 

Physically, this is a result of two compet- 
ing processes. One is the amplification pro- 
cess due to the interaction between the atom 
and the cavity. The other is the dissipation 
process due to the field-environment interac- 
tion. These two processes compete against 



each other as the system evolves. At first, 
the amplification process dominates the dis- 
sipation process and the entanglement of the 
system increases with time. Later, the dissi- 
pation process is in control and the entan- 
glement decreases with time. Finally, the 
atom- field system is disentangled as one can 
see from Fig.3. One can also use the nega- 
tivity to evaluate the degree of entanglement 
of the two-qubit system ( e.g., see Ref.[24j). 
For a 2 x 2 quantum system the concurrence 
and negativity are both good entanglement 
measures. In fact it has been proved that 
th e negativity of a state is always lager than 
y/(l -C) 2 + C 2 - (1 - C) and smaller than 
C, where C is the concurrence of a state 23 . 
For a 2 x 3 quantum system the concurrence 
is no t ap plicable while the negativity is still 
valid m 



IV. PURITY OF SYSTEM 

In this section, we investigate the purity 
of the system by employing the linear en- 
tropy. Many protocols in quantum informa- 
tion processing require pure, maximally en- 
tangled quantum states. For example, quan- 
tum teleportation often relies heavily on the 
purity and entanglement of the initial state. 
However, an pure and entangled quantum 
system usually becomes mixed and/or less 
entangled under the influence of decoherence. 
Here, we adopt the linear entropy to quantify 
the mixedness of a state defined by 



S(p) = l-Tr(p 2 ). 



(25) 



Generally, if p is the density matrix of a pure 
state, s = 0, otherwise s > 0. It has also 
been proved that a bipartite mixed states is 
useless for quantum teleportation if its linear 
entropy exceeds 1/2 for a two qubits s yste m. 
The purity of the atom-field system is[25| 



S(p) = 1 - Tr(p 2 ) = 1 - Tr F {Tr A (p 2 )} 
= 1 -Tr F {(p 00 + p n ) 2 } 

= ^(i-\U(t)\ 2 ), 



(26) 



2 t 

i 




1 



FIG. 2: The concurrence is plotted MS c\ function 
of g and t with a = 1. Upper panel: The decay 
rate constant k is zero. Lowe panel: The decay 
constant k is 0.0.5. 

where we have used the fact 

Tr F (\a + (t))(a_(t)\) = (a_(t)\a + (t)). (27) 

In order to analyze the purity of the two- 
level atom at time t, we first trace out the 
field's variables and obtain 

PA{t) = \{l + Re[f + {t)\{t)\}\V{U 

+l{l-Re[f + (t)\(t)]}\l)(l\ (28) 

+ iIm[f f m (\m\-\i)(U), 



FIG. 3: The concurrence is plotted function 
of t with g = 1 for k = 0(solid line), k = 0.1 (dash 
line) and k = l(dot line). 

with Re and Im denoting the real and the im- 
age part of a complex number, respectively. 
The linear entropy of the atom is 

s A (PA) = lii-\Mtmn (29) 

Finally, we investigate the purity of the 
field. After tracing out the atom's variables, 
the density matrix of the field is 

p F (t) = Poo(0+Pn(0 (30) 
= ^[\ a+ (t))(a + (t)\ + \a4t))(a4t)\]. 

One can write the linear entropy of the field 
as follows 

S F (p F )= 1 -[l-\\(t)\ 2 ]. (31) 

By using the explicit expression of the func- 
tions |/ + (t)| 2 and |A(t)| 2 , we find that the lin- 
ear entropies S, Sa, and S F are independent 
of the parameter a. The linear entropy is 
plotted as a function of time t for several val- 
ues of k/g in Figs. (4-6). From Fig. 4, we can 
see that the linear entropies S, Sa, and S F in- 
crease with time and then reach the maximal 
value 0.5. However, the situation is different 
when the parameter k/g increases. In Fig. 5 
and Fig. 6, it is easy to see that the maximal 
values of S F (dotted line) is less than 0.5. In 
general, the linear entropy Sa is larger than 
the linear entropy S. As the parameter k/g 
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FIG. 4: The linear entropy is plotted as a FIG. 5: The linear entropy is plotted as a 
function of t with k/g = 0.3 for S^solid line), function of t with k/g = 1 for S^solid line), 



(dotted line) and ^(dashed line). 



Sf (dotted line) and Sa (dashed line). 



increases, the linear entropy Sp decreases sig- 
nificantly(see the dotted line in Fig.6). This 
result is similar to that obtained in Ref . [2o| . 
The atomic and field coherence loss can be 
measured by the linear entropy Sa and Sf 
25[. Comparing Figs. (4-6), one can see that 



the atomic coherence Sa and the field coher- 
ence Sp decrease with the increase of the de- 
cay rate k. However, the atom-field system 
preserves its linear entropy S as the system 
evolves. For example, the linear entropy of 
the atom-field system is 0.5 if t > 5(see the 
solid lines in Figs. (4-6)) which is a result of 
two competing processes. One is the pro- 
cess due to the interaction between the atom 
and the field. The other is the process due 
to the field-environment interaction. As we 
have pointed out, these two processes com- 
pete against each other during the whole evo- 
lution. The linear entropy S first increases 
with time t and then reaches a plateau. We 
note that the similar features have also been 
observed in Ref.[l8| where a Cooper-pair box 
is put into a phase-damped cavity. 



V. CONCLUSIONS 

In the present paper, we have investigated 
entanglement dynamics of a quantum system 
consisting of one two-level atom coupled to a 
dissipative cavity with a strong classical driv- 




t 



FIG. 6: The linear entropy is plotted as a 
function of t with k/g = 6 for S^solid line), 
Sf (dotted line) and Sa (dashed line). 

ing field. An analytical expression of the time 
evolution density matrix operator for the sys- 
tem is found and used to study the entangle- 
ment and the purity of the system. There are 
two competing processes. One is the amplifi- 
cation process due to the interaction between 
the atom and the cavity. The other is the dis- 
sipation process due to the field-environment 
interaction. Our calculation shows that there 
is an optimal interaction time to maximize 
the entanglement of the atom- field system in 
the presence of dissipation when the coupling 
constants and the frequencies of the two-level 
atom, the cavity, and the classical driving 
field are given. The approach adopted in the 
present paper may be extended to systems 
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of two or more two-level atoms in dissipative 
cavities. 
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